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Abstract. One of the observational evidences in support of the thermonuclear runaway model for 
Q . the classical nova outburst relies on the accompanying nucleosynthesis. In this paper, we stress 

^0 ' the relevant role played by nucleosynthesis in our understanding of the nova phenomenon by 

— ^ \ constraining models through a comparison with both the atomic abundance determinations from 

1*"^ ■ the ejecta and the isotopic ratios measured in presolar grains of a likely nova origin. Furthermore, 

the endpoint of nova nucleosynthesis provides hints for the understanding of the mixing process 
, responsible for the enhanced metallicities found in the ejecta, and reveals also information on the 

^ ■ properties of the underlying white dwarf (mass, luminosity...). 

\ We discuss first the interplay between nova outbursts and the Galactic chemical abundances: 

OO . Classical nova outbursts are expected to be the major source of ^^C, '^N and '^O in the Galaxy, 

^""l ' and to contribute to the abundances of other species with A < 40, such as ^Li or ^^Al. We describe 

the main nucleai" path during the course of the explosion, with special emphasis on the synthesis 
^—^ . of radioactive species, of particular interest for the gamma-ray output predicted from novae (^Li, 

■ 22Na, 26 Al). An overview of the recent discovery of presolar nova candidate grains, as well as 
, a discussion of the role played by nuclear uncertainties associated with key reactions of the NeNa- 

■ MgAl and Si-Ca regions, are also given. 

O ■ 

■ 

; GALACTIC ALCHEMY: THE INTERPLAY BETWEEN NOVA 

; OUTBURSTS AND THE GALACTIC ABUNDANCES 

. 5^ ! The high peak temperatures achieved during nova explosions, T^^^^ ~ (2 — 3) x 10^ K, 

suggest that abundance levels of the intermediate-mass elements in the ejecta must be 
significantly enhanced, as confirmed by spectroscopic determinations in well-observed 
nova shells. This raises the issue of the potential contribution of novae to the Galac- 
tic abundances, which can be roughly estimated as the product of the Galactic nova 
rate, the average ejected mass per nova outburst, and the Galaxy's lifetime. This order 
of magnitude estimate points out that novae scarcely contribute to the Galaxy's overall 
metallicity (as compared with other major sources, such as supernova explosions), nev- 
ertheless they can substantially contribute to the synthesis of some largely overproduced 
species (see Table 1, for a sample of publications addressing nucleosynthesis in classical 
novae). Hence, classical novae are likely sites for the synthesis of most of the Galactic 
^^C, ^^N and ^^O, whereas they can partially contribute to the Galactic abundances of 
other species with A < 40, such as "^Li, ^^F, or ^^Al [83, 39]. 

Overproduction factors, relative to solar, corresponding to hydrodynamic calculations 



TABLE 1. A sample of publications in refereed journals addressing nucleosynthesis in classical 
novae 



Reference 


Model category 


Range of nuclei 


Amould & N0rgaard (1975) [4] 


Parametric, 1 zone 


^He, ^Li, B, C, N, 


Arnouldetal. (1980) [5] 


Parametric, 1 zone 


H-Ar 


Boffin etal. (1993) [8] 


Parametric, 1 & 2 zones 


^He, ^Be, ^Li, ^B, 


Cocetal. (1995) [13] 


Semianalytic, 1 zone 


H-K 


Coc et al. (2000) [14] 


Hydrodynamic, ID 


0,F 


Glasner etal. (1997) [26] 


Hydrodynamic, 2D 


H-F 


Hemanz et al. (1996) [28] 


Hydrodynamic, ID 


^He, ■'Be, '^U 


Hernanz et al. (1999) [29] 


Hydrodynamic, ID 


"N, i^F 


Hillebrandt & Thielemann (1982) [32] 


Parametric, 1 zone 


H-Ar 


lUadis et al. (1999) [33] 


Parametric, 1 zone 


H-Ca 


Jose et al. (1997) [37] 


Hydrodynamic, ID 


Na, Mg, Al 


Jose & Hernanz (1997) [38] 


Hydrodynamic, ID 


H-Ca 


Jose & Hernanz (1998) [39] 


Hydrodynamic, ID 


H-Ca 


Jose etal. (1999) [40] 


Hydrodynamic, ID 


Ne-Na, Mg-Al 


Jose etal. (2001a) [41] 


Hydrodynamic, ID 


"N, i^F, ■'Be, 22Na, ^^Al 


Jose et al. (2001b) [42] 


Hydrodynamic, ID 


C,N,0,Al,Mg,Si,Ne 


Jose etal. (2001c) [43] 


Hydrodynamic, ID 


Si-Ca 


Kercek et al. (1998) [47] 


Hydrodynamic, 2D 


H-F 


Kerceket al. (1999) [48] 


Hydrodynamic, 3D 


H-F 


Kolb & Politano (1997) [49] 


Hydrodynamic, ID 


22Na, 26 Al 


Kovetz & Prialnik (1985) [50] 


Hydrodynamic, ID 


H-O 


Kovetz & Prialnik (1997) [51] 


Hydrodynamic, ID 


H-O 


Kudryashov & Tutukov (1995) [52] 


Parametric, 1 zone 


H-Ar 


Kudryashov et al. (2000) [53] 


Parametric, 1 zone 


H-Ar 


Lazareffetal. (1979) [54] 


Parametric, 2 zones 


C, N, O, F, Ne 


Nofaretal. (1991)[56] 


Parametric, 1 zone 


H-Al 


Politano etal. (1995) [59] 


Hydrodynamic, ID 


H-Ca 


PriaMik etal. (1978) [60] 


Hydrodynamic, ID 


H-O 


Prialnik etal. (1979) [61] 


Hydrodynamic, ID 


H-O 


Prialnik (1986) [62] 


Hydrodynamic, ID 


H-O 


Prialnik & Shara (1986) [63] 


Hydrodynamic, ID 


H-Ne 


Prialnik & Shara (1995) [64] 


Hydrodynamic, ID 


H-P 


Shara & Prialnik (1994) [69] 


Hydrodynamic, ID 


H-Mg 


Sparks etal. (1978) [73] 


Hydrodynamic, ID 


H-O 


Starrfield etal. (1972) [74] 


Hydrodynamic, ID 


H-O 


Starrfield et al. (1974a) [75] 


Hydrodynamic, ID 


H-O 


Starrfield etal. (1974b) [76] 


Hydrodynamic, ID 


H-O 


Starrfield etal. (1978a) [77] 


Hydrodynamic, ID 


^He, 'Li, C, N, O 


Starrfieldetal. (1978b) [78] 


Hydrodynamic, ID 


'Li, C, N, O 


Starrfield etal. (1992) [80] 


Hydrodynamic, ID 


H-Ar 


Starrfield etal. (1993) [81] 


Hydrodynamic, ID 


H-Ar 


Starrfieldetal. (1998) [83] 


Hydrodynamic, ID 


H-Ar 


Starrfield et al. (2000) [84] 


Hydrodynamic, ID 


H-Ar 


Starrfield etal. (2001) [85] 


Hydrodynamic, ID 


H-S 


Vangioni-Flam et al. (1980) [87] 


Parametric, 1 zone 


Ne, Al 


Wallace & Woosley (1981) [89] 


Parametric, 1 zone 


H-Al 


Wanajo et al. (1999) [90] 


Semianalytic, 1 zone 


H-Ca 


Weiss & Truran (1990) [91] 


Parametric, 1 zone 


H-Ca 


Wiescher et al. (1986) [92] 


Parametric, 1 zone 


H-Ar 
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FIGURE 1. (Left) Mean overproduction factors, relative to solar, in the ejecta of a 1.15 Mq CO novae. 
(Right) Same for a 1.35 Mq ONe novae. Figures are based on hydrodynamic calculations reported in [39]. 



of nova outbursts on top of a 1.15 Mq CO and a 1.35 Mq ONe white dwarf, are 
shown in Figure 1. Because of the lower peak temperatures achieved in CO models, 
and also because of the lack of significant amounts of seed nuclei in the NeNa-MgAl 
region, the main nuclear activity in CO novae does not extend much beyond oxygen, 
as seen from the overproduction plot. In contrast, ONe models show a much larger 
nuclear activity, extending up to silicon (1.15 Mq ONe) or argon (1.35 Mq, ONe). 
Hence, the presence of significantly large amounts of intermediate-mass nuclei in the 
spectra, such as phosphorus, sulfur, chlorine or argon, may reveal the presence of an 
underlying massive ONe white dwarf. Another trend derived from the analysis of the 
nucleosynthesis accompanying nova outbursts is the fact that the 0/N and C/N ratios 
decrease as the mass of the white dwarf (and hence, the peak temperature attained during 
the explosion) increases. 



Abundance Determinations in the Ejecta from Novae 

In order to constraint the models, several works have focused on a direct comparison 
of the atomic abundances inferred from observations of the ejecta with the theoretical 
nucleosynthetic output (see [39, 83], and references therein). Despite of the problems 
associated with the modeling of the explosion [86], such as the unknown mechanism 
responsible for the mixing between the accreted envelope and the outermost shells of the 
underlying white dwarf [9, 21], or the difficulties to eject as much material as inferred 
from observations [71], there is an excellent agreement between theory and observations 



as regards nucleosynthesis (i.e., including atomic abundances -H, He, C, O, Ne, Na-Fe-, 
and a plausible endpoint for nova nucleosynthesis). In some cases, such as for PW Vul 
1984, the agreement between observations and theoretical predictions (see [39], Table 
5, for details) is really overwhelming. The reader is referred to [24] for an extended list 
of abundance determinations in the ejecta from novae, and to [68, 88] for recent efforts 
to improve the abundance pattern for QU Vul 1984 and V1974 Cyg 1992, respectively. 

Since the nuclear path is very sensitive to details of the evolution (chemical compo- 
sition, extend of convective mixing, thermal history of the envelope...), the agreement 
between inferred abundances and theoretical yields not only validates the thermonuclear 
runaway model, but also poses limits on the (yet unknown) mixing mechanism itself: for 
instance, if mixing occurs very late in the course of the explosion, the accumulation of 
larger amounts of matter in the envelope will be favored (since the injection of signifi- 
cant amounts of the triggering nucleus ^^C will be delayed). Hence, one would expect to 
end up with a more violent outburst, characterized by a higher T^^^^, exceeding in some 

cases 4 x 10^ K, and, as a result, a significant enrichment in heavier species, beyond 
calcium, in the ejecta from novae involving very massive white dwarfs, a pattern never 
observed so far. 



Presolar Grains: Gifts from Heaven 

Infrared [22, 24] and ultraviolet observations [70] of the temporal evolution of nova 
light curves suggest that novae form grains in the expanding nova shells. Both CO 
and ONe novae behave similarly in the infrared right after the outburst. However, as 
the ejected envelope expands and becomes optically thin, such behavior dramatically 
changes: CO novae are typically followed by a phase of dust formation corresponding 
to a decline in visual light, together with a simultaneous rise in the infrared emission 
[23, 25]. In contrast, it has been argued that ONe novae (that involve more massive 
white dwarfs than CO novae) are not so prolific producers of dust as a result of the lower 
mass, high-velocity ejecta, where the typical densities can be low enough to enable 
the condensation of appreciable amounts of dust. Hints on the condensation of dust 
containing silicates, silicon carbide, carbon and hydrocarbons have been reported from 
a number of novae (see [24] for a recent review). 

Up to now, the identification of presolar nova grains, presumably condensed in the 
shells ejected during the explosion, relied only on low ^'^Ne/^^Ne ratios (attributed to 
^^Na decay), but quite recently five silicon carbide and two graphite grains that exhibit 
isotopic signatures characteristic of nova nucleosynthesis have been identified [1, 2]. 
They are characterized by very low ^^C/^^C and ^"^N/^^N ratios, •^'^Si excesses and close- 
to- or slightly lower-than-solar ^^Si/^^Si ratios, high ^^Al/^^Al ratios (determined only 
for two grains) and low ^°Ne/^^Ne ratios (only measured in the graphite grain KFBla- 
161). Such a promising discovery provides a much valuable source of constraint for nova 
nucleosynthesis (since contrary to the atomic abundance determinations derived from 
nova ejecta, measurements provide more accurate isotopic ratios) and opens interesting 
possibilities for the future. 

Theoretical isotopic ratios for a variety of nuclear species, ranging from C to Si, based 



on hydrodynamic computations of the nova outburst, have been reported by different 
authors [82, 42, 44, 45]. A more detailed analysis, which focuses on the different 
chemical pattern expected for CO and ONe novae, will be presented elsewhere [46]. 



SYNTHESIS OF RADIOACTIVE NUCLEI DURING NOVA 

OUTBURSTS 

Among the isotopes synthesized during classical nova outbursts, several radioactive 
species deserve a particular attention. Short-lived nuclei, such as ^^'^^O and ^^F (and 
to some extent ^^N) have been identified as the key isotopes that power the expansion 
and further ejection during a nova outburst through a sudden release of energy, a few 
minutes after peak temperature [74]. Other isotopes have been extensively investigated 
in connection with the theoretical gamma-ray output from novae [11, 12, 55]. Hence, ^^N 
and ^^F are responsible for the predicted prompt y-ray emission [29] at and below 511 
keV, whereas ^Be and ^^Na [27], which decay much later, when the envelope is optically 
thin, are the sources that power line emission at 478 and 1275 keV, respectively. ^^Al 
is another important radioactive isotope that can be synthesized during nova outbursts, 
although only its cumulative emission can be observed because of its slow decay. 

We will briefly focus on the corresponding nuclear paths leading to the synthesis of 
the abovementioned gamma-ray emitters, with special emphasis on the nuclear uncer- 
tainties associated with the relevant reaction rates. We refer the reader to [30] for a 
review of the current theoretical predictions of the gamma-ray output from novae and 
of the chances for a nearby future detection using spacecrafts such as INTEGRAL (see 
also [31]). As for a comprehensive summary of the main uncertainties affecting nuclear 
reaction rates for nova temperatures, the reader is referred to [15]. Other recent attempts 
to fully analyze the impact of nuclear physics uncertainties in nova nucleosynthesis in- 
volve parametric one-zone calculations with temperature and density profiles extracted 
from hydrodynamic models [34, 35], as well as Monte Carlo simulations [72]. We would 
like to stress, however, that results based on parametric calculations usually tend to over- 
estimate the effect of a given reaction rate uncertainty, as compared with the outcome 
from hydrodynamic tests. Hence, although these simplified techniques provide a very 
valuable tool to potentially identify key reactions through an extraordinary large number 
of tests, final conclusions have to be confirmed through hydrodynamic calculations. 

^Be-^Li 

■^He(a,7)^Be is the main reaction leading to ^Be synthesis (with ^He coming from 
both the accreted amount plus the contribution in place from ^H(p,e+ Ve)^H(p,7)^He), 
which is transformed into ^Li by means of an electron capture, emitting a y-ray photon 
of478keV [27]. 

Its production in classical nova outbursts has been very controversial. Results from 
the first pioneering calculation, in the framework of a simple parametric model [4], 
were confirmed by early hydrodynamic simulations [77], assuming however an enve- 



lope in-place (the first hydrodynamic nova models that properly included the onset of 
the accretion phase were not available until the 80s), thus neglecting the impact of the 
accretion phase on the evolution. These results were refuted later on, in terms of para- 
metric one/two zone models [8], pointing out the critical role played by the ^C(p,7) 
reaction, not included in all previous works (i.e., [4, 77]), and claiming therefore for 
an unlikely synthesis of ^Li in novae. The scenario was recently revisited by [28, 39], 
who performed new hydrodynamic calculations, taking into account both the accretion 
and explosion stages, and a full reaction network (including ^C(p,7)). These studies con- 
firmed that the Be-transport mechanism [10] is able to produce a large overproduction 
of ^Li in nova explosions. 

Among the issues that affect ^Li synthesis in novae, one the most critical ones is the 
final amount of ^He that survives the early rise in temperature when the thermonuclear 
runaway ensues. In particular, the different timescales to reach T^^^^ achieved for CO 

and ONe novae, which deeply depend on the initial ^^C content in the envelope, lead 
to a larger amount of ^Be in CO novae (which survives destruction through ^Be(p,7)^B 
because of the very efficient inverse photodisintegration reaction on ^B). No relevant nu- 
clear uncertainties in the domain of nova temperatures affect the corresponding reaction 
rates. 

The puzzling ^Li synthesis in novae, suggested for many years but elusive up to 
now, has been apparently confirmed for the first time in the spectra of V382 Vel (Nova 
Velorum 1999) [18], for which an observed feature compatible with the doublet at 6708 
A of Lil has been reported. Although this identification is not yet confirmed (there is 
no unambiguous interpretation for such a feature), the lithium abundances inferred are 
fully compatible with the theoretical upper limits given for a fast nova (i.e., [28, 39]). 

The potential contribution of classical novae to the Galactic ^Li content turns out to 
be rather small (i.e., less than 15%), even if the ^Li yield from the most favorable case, 
a massive CO nova [28], is considered. Nevertheless, a nova contribution is required to 
match the ^Li content in realistic calculations of Galactic chemical evolution [66, 67]. 



The potential role of Na for diagnosis of nova outbursts was first suggested by [1 1]. 
It decays to a short-lived excited state of ^^Ne (with a lifetime of t = 3.75 yr), which 
de-excites to its ground state by emitting a y-ray photon of 1.275 MeV. Through this 
mechanism, nearby ONe novae within a few kiloparsecs from the Sun may provide 
detectable /-ray fluxes [30]. Several experimental verifications of this 7-ray emission 
at 1.275 MeV from nearby novae have been attempted in the last twenty years, using 
balloon-borne experiments and detectors on-board satellites such as HEAO-3, SMM, or 
CGRO, from which upper limits on the ejected ^^Na have been derived. In particular, the 
observations performed with the COMPTEL experiment on-board CGRO of five recent 
Ne-type novae [36], as well as observations of standard CO novae, have led to an upper 
limit of 3.7 X 10"^ M,,for the 22Na mass ejected by any nova in the Galactic disk. A 
limit that poses some constraints on pre-existing theoretical models of classical nova 
explosions. 



Synthesis of Na in novae, extensively investigated in the last two decades [32, 92, 
94, 91, 58, 59, 13, 52, 49, 39, 83, 40, 90, 53, 84, 41, 85], proceeds through different 
reaction paths. In the ^'^Ne-enriched envelopes of ONe novae [40], it takes place through 
^''Ne(p,7)^^Na, followed either by another proton capture and a j8+-decay into ^^Na 
(i.e., ^^Na(p,7)^^Mg(/3+)^^Na), or decaying first into ^^Ne before another proton capture 
ensues (i.e., ^^Na(j8+)^^Ne(p,7)^^Na). Other potential channels, such as proton captures 
on the seed nucleus ^^Na, play only a marginal role on ^^Na synthesis because of the 
much higher initial ^^Ne content in such ONe models. As for the main destruction 
channel at nova temperatures, ^^Na(p,7)^-^Mg competes favorably with ^^Na(/3+)^^Ne. 
Nuclear uncertainties strongly affect the rates for ^^Na(p,7)^^Mg and ^^Na(p,7)^^Mg 
[40], which translate into an uncertainty in the final ^^Na yields (and ultimately on the 
maximum detectability distance of the 1275 keV line expected from nova outbursts). 
Advances to reduce the uncertainty affecting the ^^Na(p,7) rate have been recently 
achieved with the DRAGON recoil separator facility at TRIUMF [57, 7]. 

26 Al 

Several isotopes should be considered as potential seeds for ^^Al synthesis: in par- 
ticular, ^"^'^^Mg and to some extent ^^Na and ^^'^^Ne [37]. The main nuclear reac- 
tion path leading to ^^Al synthesis, also investigated in a large number of papers 
[54, 5, 32, 17, 92, 94, 93, 91, 56, 58, 69, 59, 13, 52, 49, 39, 83, 40, 90, 53, 84, 41, 85], 
is given by ^^Mg(p,7)^^Al(j8+)^^Mg(p,7)^^AF, whereas destruction is dominated by 
26Al«(p,7)26si. 

A significant nuclear uncertainty affects the ^^Al(p,7)^^Si rate [40], which translates 
into an uncertainty in the expected contribution of novae to the Galactic ^^Al content. 
A critical issue in order to estimate this contribution is the initial composition of the 
ONe white dwarf. Whereas calculations by Starrfield et al. [79, 83, 84, 85] assume a 
core composition based on hydrostatic models of carbon burning nucleosynthesis [3], 
rather enriched in ^^Mg (with a ratio ^^0:^*^Ne:^^Mg around 1.5:2.5:1), we adopt more 
recent values taken from stellar evolution calculations of intermediate-mass stars [65], 
for which the ^"^Mg content is much smaller (^^0:^°Ne:^^Mg is 10:6:1). Calculations 
based on the new ONe white dwarf composition [37, 39, 40, 41], suggest that the 
contribution of novae to the Galactic ^^Al abundance is rather small (i.e., less than 15%), 
in good agreement with the results derived from the COMPTEL map of the 1809 keV 
26 Al emission in the Galaxy (see [20]), which points towards young progenitors (type II 
supemovae and Wolf -Ray et stars). 

18p 

The predicted gamma-ray emission from novae at and below 511 keV at early 
epochs after the explosion is basically driven by the amount of ^^F present in the en- 
velope [55, 27, 29]. The synthesis of ^^F is powered by ^^0(p,y)^^F, followed either 
by i'^F(p,7)i^Ne(j8+)i^F or by i'^F(j8+)i'^0(p,7)i^R The dominant destruction chan- 



nel is ^^F(p,a)^^0 plus a minor contribution from ^^F(p,7)^^Ne. The effect of the nu- 
clear uncertainties associated with some of the rates (i.e., ^^F(p,ct)^^0, ^^F(p,7)'^Ne, 
i'^0(p,a)i4Nand ^^0{p,r)^^¥) [29, 14, 15] is, in this case, quite remarkable: they trans- 
late into a large uncertainty in the expected ^^F yield, and therefore, in the corresponding 
gamma-ray flux and maximum detectability distance. Advances to reduce this uncer- 
tainty have been reviewed during the Classical Nova Conference [6, 19] and elsewhere 
[16], and involve several nuclear physics experiments performed in Oak Ridge (USA) 
and Orsay (France). 



The endpoint of nova nucleosynthesis 

In agreement with the chemical pattern derived from detailed observations of the 
ejecta, the theoretical endpoint for nova nucleosynthesis is limited to A < 40 (i.e., 
calcium), in agreement with current theoretical nucleosynthetic estimates, provided 
that the temperatures attained in the envelope during the explosion remain limited to 
V^-(2-3)x108k. 

The nuclear activity in the Si-Ca region has been scarcely analyzed in detail in the 
context of classical nova outbursts [83, 33, 90, 43]. It is powered by a leakage from 
the NeNa-MgAl region, where the activity is confined during the early stages of the 
explosion. The main reaction that drives the nuclear activity towards heavier species 
(i.e., beyond S) is mainly 30p(p,7)3is, either followed by ^^S(p,Y)^^Cl(j5+)^^S, or by 
^^S(j8+)^^P(p,7)-^^S [43]. The ^°P(p,7) rate is based only on Hauser-Feshbach estimates, 
which can be rather uncertain at the domain of nova temperatures. To test the effect of 
this uncertainty on the predicted yields, we have performed a series of hydrodynamic 
calculations [43], modifying arbitrarily the nominal rate. Hence, for a high ■''^P(p,7) rate 
(i.e., 100 times the nominal one), the final -^"Si yields are dramatically reduced by a 
factor of 30, whereas for a low ^*^P(p,7) rate (i.e., 0.01 times the nominal one), the final 
3*^Si yields are slightly increased by a factor of 5, whereas isotopes above silicon are 
reduced by a factor of ~ 10, with dramatic impact on theoretical estimates ionvolving 
both the composition of the ejecta and of presolar grains. 

The impact of nuclear reaction rate uncertainties on nucleosynthesis calculations 
points out the need of accurate nuclear physics inputs, and stresses the role played by 
classical novae as perfect laboratories for nuclear astrophysics. 
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